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Abstract 

We study the relation between the coefficients of Taylor series and Kapteyn 
series representing the same function. We compute explicit formulas for 
expressing one in terms of the other and give examples to illustrate our 
method. 
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1. Introduction 



Series of the form 26 



^a^Jn+u [{n + v)z\, (1) 



n=0 



and 



J^^'^^J^+n + + 2n) z] J,+„ + P + In) z] , (2) 

ra=0 

where E C and Jn{') is the Bessel function of the ffist kind, are called 
Kapteyn series of the first kind and Kapteyn series of the second kind re- 
spectively. 
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Kapteyn series have a long history, going back to Joseph Louis de La- 
grange's 1771 paper Sur le Probleme de Kepler IJ], where he solved Kepler's 
equation [sj 

M = E-esm{E), (3) 

using his method for solving implicit equations [13] (now called Lagrange 
inversion theorem) and obtained ^ 

Jn—l 



E{M) = M + y — —— sin" (M) . 



n=l 



Here M is the mean anomaly (a parameterization of time) and E is the 
eccentric anomaly (an angular parameter) of a body orbiting on an ellipse 
with eccentricity e. 

In 1819 Friedrich Wilhelm Bessel published his paper Analytische Aufldsung 
der Kepler' schen Aufgabe flJ], where he approached using a different 
method. First of all he observed that the function g{M) = E{M) — M defined 
imphcitly by = esin {g + M) is 27r— periodic and satisfies 5^(0) = = g{7i). 
Hence, g{M) can be expanded in a Fourier sine series 

oo 

g{M) = Y^bn sin (nM) , 



n=l 



where 



g{M) sm{nM)dM 
cos{nM) 



g{M)- 



n 



2 f 

H / cos{nM)dg 

™i 



TT 

— [ cos(nM)d(E - M) 
Tin J 



/ COS [uiE-e sin E)] dE - — I 
J rcnj 



2 

nn 



cos{nM)dM 



2 



and hence 

, 2 

bn = 

irn 



J cos {nE — ne sin E) dE. 



He then introduced the function Jn{z) defined by 

Jniz) = - J COS (nE - z sin E) dE, neZ (4) 



which now bears his name and obtained 



2 

E{M) = M + nine) sin (nM) . (5) 



'n 

71=1 



Bessel's work on (j4]) was continued by other researchers including Lommel, 



who defined the Bessel function of the first kind by 18 

- (-1)" /£ 
'n\T (u + n + l) \2 



J^W = E ..,rLLn y ■ "'^C, (6) 



n=0 



where r(-) is the Gamma function. 

In 1817, Francesco Carhni jij found an expression for the true anomaly 
V (an angular parameter), defined in terms of E and e by 



/v\ 1 + e E 
tan - = A / tan — 

Carlini's expression reads 0] 

oo 

v = M + sin (nM) , 

where 

2 °° 



n 

m=0 



with e = ■ The problem considered by Carlini was to determine the 
asymptotic behavior of the coefficients Bn for large values of n [9|. The as- 
tronomer Johann Franz Encke drew Carl Gustav Jacob Jacobi's attention to 



3 



the work of Carlini. In 1849, Jacobi published a paper improving and cor- 
recting Carhni's article [lo| and in 1850 Jacobi published a translation from 
Italian into German [sl, with critical comments and extensions of Carlini's 
investigation. 

Bessel's research on series of the type (j5]) was continued by Ernst Meissel 



24j | in his papers 19|], 20| and by Willem Kapteyn (not to be confused with 



his brother Jacobus Cornelius Kapteyn 23|) in the articles ll| and 12 
Most of the early work on Kapteyn series, together with their own results, 
can be found in the books by Niels Nielsen |22|, Chapter XXII] and George 



Neville Watson j26|. Chapter 17]. 

In recent years, there has been a renewed interest on Kapteyn series, par- 
ticularly from researchers in the fields of Astrophysics and Electrodynamics 
(see 17| for a review of current applications). In ItlSj], Ian Lerche and Robert 



C. Tautz studied the Kapteyn series of the second kind 



Si{a) = y^n' 



,4 t2 



J^ (ria) 



n=l 



and derived the formula 



S'i(a) 



(64 + 592a2 + 4720^* + 27a^ 



256(1 - a' 



They continued their investigations in p^, where they outlined a way for 
calculating more general Kapteyn series of the form 



5*1 (m, a) = 



na] 



m 



0,1, 



(7) 



The purpose of this paper is to describe a method for computing the 
coefficients in the Taylor series of functions defined by Kapteyn series of the 
first (II]) and second ([2]) kind. As an example, we will show a closed-form 
expression of ([7]) valid for all values of m. 



2. Kapteyn series of the first kind 

We begin by considering functions expressed as Kapteyn series of the first 
kind. 



4 



Theorem 1. Suppose that 



m=0 



and 



7 (z) = ^a^J^+n [{t^ + n) z] 



(8) 



(9) 



n=0 



where both series converge absolutely for z in some domain Q. Then, we have 

(10) 



LfJ 

y^ys,mb 



s—2m 



m=0 



and 



for all values of with 



LfJ 

E 

m=0 



-1^ 



BTiv + n-k + l] 



u + n — 2k 



and 



Vn,k 



l{u + n- 2kf T{v + n-k) [ 2 



A;! 



p + n 



u+n 



v+n-2k+l 



Proof. Let z eVL.To prove ([TT]), we use ^ and ([9]), to get 

oo oo 

= z'^f (z) = J]<J.+„ [{ly + n) z] 



m=0 

and from (El) we obtain 



n=0 



oo oo 



•-11 



n=0 m=0 



m!r (i/ + n + m + 1) 



(zy + n) 2 



v+n+2m 



or 



oo oo 



-1^ 



s=0 



n=0 m=0 



m\V + n + m + 1) \ 2 



z/ + n 



u+n+2'm 



Z 



n+2m 



(11) 



:i2) 
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Setting n + 2m = s on the right-hand side of (fT^ , we have 

(-1) (^s-2m (v + s-2m^ 



^ ^ ^m!r (z/ + s - m + 1) 

s=0 s=0 m=0 ^ ' 

from which (fTTj) follows. 
We have 

(u + 2k) fu + p + s 



(s — fc) (z^ + s + /c) \ p — k 
and hence 

^^i2ij-fcW2s,s-i = 0, k S. 
j=k 

When k = s, we get 

1 fiy + 2sY^'' 

- r(^ + 2. + l) [~^) 

u+2s 



X (z/ + 2.) r (z. + 2.) (^-^J =1 

and therefore 

s 

'y^^U2jJ-kV2s,s-j = ^k,s- 
j=k 

Thus, 

s s s j 

''^^V2s,mb2{s-m) = ^^t^2s,s-j^2j = ^^^2s,s-j ^^^2j,ma2(j-m) 
m=0 J=0 J=0 m=0 

s j s s 

= ^V2s,s-j^U2j,j-ka2k = ^Ci2k'^^2j,j-kV2s,s-j = 
j=0 k=0 k=0 j=k 

s 

^^^2fc'^fc,s = (^2s 
k=0 

and a similar computation holds for 62^+1 and a2fc+ii proving flTU]) . 
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Remark 2. Formula ^^) appeared in I2d . 17.5 (6)], hut it contains a small 
mistake because there is a factor of | missing in the denominator. The result 
was also published in 7.10.2 (29)], but there is also a misprint there, since 
the factor hn-2m is missing in the sum. 

In f^], we computed formulas for hg and a'^ when z/ = 0. 

In order to find the coefficients hg for a particular choice of = and aj^, 
we need the foUowing result. 



Lemma 3. Let r,m & N. Then, we have 

r 



k=0 

Proof. We have 



dt"-- 



t=o 



sinh"(t) 



e — e ^ 



2'^^( A; 

fc=0 ^ 



k (r-2k)t 



-Ife 



Since 



we obtain, 



"m „at 



dt'' 



d"" sinh'^(t) 



dr' 



k=0 ^ ^ 



(13) 



Example 4. LeVs consider the special case u = and a'^ = v?^ . Then, ( flil) 
gives 

^^(^) = 7^E(-l)'(l)(^-2^)"''- (14) 

k=0 ^ ^ 

Since the terms in the sum are symmetric with respect to k = [|J , we can 
write 

hs{p) ^ ^' 



k=0 ^ ' 



s+2p 



unless s = = p, in which case we have 

bo{0) = 1. 
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Thus, we have 



k=0 ^ ^ 



: + 2p 



with 



1, s = = p 
i, otherwise 



(15) 



From [T^) . we obtain 



, , , es,p d'+^Psmh'{t) 



or 



sinh(t) 



(16) 



where [t^]G{t) denotes the coefficient of in the McLaurin series of G{t). 



From fTS\) we get 






bsiO) = 




0, ^'.(1) 


bs{2) = 


1 
2 


(^ + 1)4- 


bs{3) = 


1 




2 


(^ + 1)6- 


bs{^) = 


1 
2 


(^ + 1)8- 



1 

2 



) - 



360 



45360 



5443200 



// we define 



s=0 



n=0 

it follows from fTd\) that 

00 



s=0 



sinh(t) 



8 



^ sinh(t) 



-{2p+l) ■ 



(17) 



Using P4 ), we g'ei 
/o (^) = 

h (^) = 

anc? in general 



2- z 
2(1 - z 
;2(9^ + 1) 
2(1-^)' 



, /3(^) 



,4 ' 



2(1-2 
2 (255^2 ^ 54^ + 1) 



2(1-^: 



10 



(11025^3 ^ 4131^2 ^ 243z + 1) 



/n(2;) 



Z Pn~l(z) 



2(l-.j 

where Po{z) = 1 and Pniz) is a polynomial of degree n for n = 1,2, 
Using the relation jM, 1 7. 33] 



we /inc? that the polynomials Pn{z) satisfy 

Pn+i{z) = P':i{z) + ^ [(1 - 6n) z' + (6n - 4)2 + 3] P'^{z) 

2+1 2 + 1 -' 

+ r9n222 + (9n + 1)2 + ll F„(2). 

2 + 1 

3. Kapteyn series of the second kind 

We now consider functions expressed as Kapteyn series of the second 
kind. 



Theorem 5. Suppose that 



(18) 



m=0 



9 



and 



(^) = + + + 2n) ^] J,+„ [(/i + /V + 2n) 



n=0 



(19) 

where both series converge absolutely for z in some domain Q. Then, we have 



2k+l 



k=0 



k=0 



(20) 



and 



for all fi, z/, with 



fc=0 



k=0 



(21) 



(/i + + 2A:)' r (/X + A; + 1) r (;/ + A; + 1) 
(// + z/ + s + A;) (^ + 1/ + 2s) 



;Lt + zy + s + A; 
s — k 



+ u + 2s 



fi+u+2k 



and 



s,k 



■1) 



s+fc 



r(^ + s + i)r(z/ + s + i)v s-k 



^ + V + 2s\ ( + 1/ + 2k 



Proof. Let z E Q. Using the formula [26|, 5.41] 



J, (z) J. (z) = J2 



k=0 



T{fi + k + l)T{u + k + l) 



fi + iy + 2k\ fzY+^+'ik 
k 



in dH]), we get 

oo oo 



-1)' 



r(/i + n + A; + l)r(z/ + n + A; + l) 



^ + z/ + 2n + 2A; 
A; 



(/i + + 2n) z 



10 



or 



oo oo ( — 1)^ 



k 

Setting n = s — k, we have 



(-1)' 



or 



+ Z^T{^^ + s + l)T{v + s + l)i^^^ ^ \ k )\ 2 ) 



11,1/ 



Comparing with ( TTSl) . we conclude that 



r(/i + s + i)r(z/ + s + i)^' 'V ^ yV 2 



r(/i + s + i)r(^/ + s + i)^^ V k 

from which (121 p follows. 
Since 



H+U+2S 



i,_n V 



/i + z/ + 2s^ ffi + i' + 2s — 2k^ ^ 



j+p ( + v + s + p\ f n + v + 2p 
s-p J\ p-j 



= (-1) 

k=j 

^ / ^ + ^ + 2jY +-+^^+^ {p-s){f, + u + 2s) 
\fi + u + 2sj {j -s){n + iy + s+j) 



11 



we see that 



while for j = s, we have 



k=j 



<OT;: = r(/x + . + i)r(z/ + . + i) 



fi + u + 2s 



fi + 1/ + 2s 



1. 



Therefore, 



s s k 

fc=0 fc=0 j=0 

j=0 k=j j=0 

The same calculation holds for 62fc+i and c^''', proving ( l20l) . ■ 
Remark 6. Nielsen Je^J, defined Kapteyn series of the second type as 



Z\ 2 



n=0 



+ n ] z 



2 



+ n ] z 



and assuming that 



oo 

. — , / 7\ra 

F{z) = Y^B^[-) , 



m=0 



he obtained 



fc=0 



+ s - 2A; r 



A; + l r 



X 



i^ + s-k-l 

k 



ji + u + 2s 

We now have all the necessary tools to compute ([7]) 



I' + s 
2 

Bs-2k- 



k + l 
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Example 7. Let 

Then, 62^+1 = and we have 

(-1) 



hsip) 



fc=0 ^ ^ 



or 



2(s+p) 



S'mce i/ie terms in the sum are symmetric with respect to k = s, we can write 

2s 



hsip) 



2{s+p) 



unless s = = p, in which case we have 

6o(0) = 1. 

Thus, we have 



s - k) 



2{s+p) 



2k) 



2{s+p) 



where Sg^p was defined in Using ^E), we get 



hsip) 



sinh^^ (t) 



t=o 



or 



hsip) 



'S,p 



i2s + 2p)\ [t^P] 



sinh (t) 



n 2s 



(22) 



13 



and therefore 



u (n^ (2s)! l {2s + 2y. s 



62^(2) 
&2.(3) 

62.(4) 



1 (2s + 4)! s (5s- 1) 

2 42 (s!)^ 90 ' 

1 (2s + 6)!s(35s2-21s + 4) 

2 43 (s!)^ 5670 ' 

1 (2s + 8)! s (5s - 2) (35s2 - 28s + 9) 



2 44 (s!)' 
Let Qp [z) he defined by 



340200 



,2s 



n=0 



s=0 



Then, since 



, , 1\ {2s + 2py. 



we get from ( f^) t/iai 
(2p)! 



£/p(^) 



4p 



(2j9)!4^ 



^sinh (t) 



2s 



4p 



1 1 



Using (2^1 . we obtain 



9o {z) = - + 



9i 



1 



z^ (1 + z^ 



sinh^ 2' 
t^ 



92 [z) 



93 (z) 



2 2^1-4^2' - • ■ (1-4^2)2 
z^l + 37z'^ + 118z^ + 27z^) 



13 
2 



(1 -4;z2)- 

^2 (1 + 217^2 + 5036z^ + 23630^6 + 22910^^ + 2250^^°) 

(1 -4^2)^ 



14 



Remark 8. Formulas for qq (z) and gi (z) were computed by George Augus- 



tus Schott in l2d J. They were reproduced by Watson in I2d . 17.6 (2)- (3)], but 



there is a typographical mistake in the equation for gi {z) , since the denomi- 

I power, instead of the | 



nator should contain a | power, instead of the ^ printed. 
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